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SUPERSONIC AERODYNAMIC INTERFERENCE EFFECTS OF 
STORE SEPARATION, PART I: 
COMPUTATIONAL ANALYSIS OF CAVITY FLOWFIELDS 


By 

Oktay Baysal* 

1. INTROOUCTION AND BRIEF REVIEW OF EXISTING WORK 
The first phase of the investigation is devoted to the study of the 
cavity fiowfields. It is aimed at developing a predictive method based on a 
rational model of the cavity flowfield to obtain sufficiently accurate ve- 
locity and pressure field determination needed to predict store motion while 
in the region of influence of the cavity. 

The simulation efforts are being conducted in two steps. The first 
step is the simulation of a two-dimensional cavity, that is, the width is 
considered infinitely long. Such an exercise is intended to be educational 
in understanding the flow behavior around the leading edge, the front face, 
the ceiling, and the rear face of the cavity, the formation of the shear 
layer, its deflection and bridging characteristic, the possible interactions 
of the shear layer, the boundary layer and shock waves. Due to tne symmetry 
of the flowfield with respect to the h-L-plane at w/2, two-dimensional 
simulations represent the flow at this particular plane where the third 
dimension effects are eliminated (Fig. 1). When sufficient learning 
experience is achieved, and necessary confidence in the preliminary two- 
dimensional code is acquired, the final step of the cavity flowfield study, 
that is, the three-dimensional cavity, will be initiated. 

Supersonic cavity flowfields, in the broadest sense, have been 

^Assistant Professor, Department of Mechanical Engineering and Mechanics, 

Old Dominion University, Norfolk, Virginia 23508. 



investigated experimentally by Stallings et al. (Refs. 1 and 2) at NASA 
Langley Research Center. There exist empirical formulae suggested by 
earlier reports (Refs. 3, 4, and 5). More recent experimental studies have 
been conducted by both U.S. Air Force and Soviet researchers (Refs. 7-9). 
Sources of interest of these studies varied from store separation to acous- 
tics and heat transfer. 

Flow Mach numbers used vary from 0.3 to 6.0. Among the measured 
quantities were shock wave pulsations and heat transfer coefficient. 
Varfolomeyev et al. (Ref. 9) observed the vortices mentioned above in their 
experiments. McGregor et al. (Ref. 6) proposed a drag coefficient relation 
for the cavity and concluded that the effect of the wind tunnel was on the 
strength of the pressure waves but not on their patterns. Wind tunnel ef- 
fects diminished with increases in flow Mach number. Sinha et al.(Ref. 7) 
concluded that the free shear region showed self-similarity and the recircu- 
lating flow had the character of a wall jet in the reverse direction super- 
imposed on a forward moving free shear layer. With the more efficient 
application of numerics and advanced computational resources, analytical 
simulations of the cavity flowfields started appearing in the literature 
(Refs. 10-16). Although these efforts are presently restricted by assump- 
tions, such as, inviscid, laminar, incompressible, two-dimensional, or only 
certain length-to-height ratios, more complete studies are in the develop- 
ment stage at present. 

2. DESCRIPTION OF PHYSICAL AND MATHEMATICAL MODELS 
Flow Development ; 


With the separation of the store, a cavity is created in the oncoming 



supersonic flow. This causes a significant increase in the fluid-dynamic, 
fluid-resonant drag, and generate fluid-elastic oscillations. The forward 
corners generate the velocity and pressure fluctuations which are amplified 
during the expansion and the convection of the flow downstream. The oncom- 
ing viscous flow will separate from the cavity forward edge, generating the 
highly vortical shear layer around the top and counter rotating eddies near 
the bottom corners as a result of the adjustment of the shear and pressure 
forces acting on the recirculating fluid to the no-slip condition imposed by 
the walls. Viscous forces acting near body surfaces can reduce the momentum 
of the fluid in a thin boundary layer region, so that when acted upon by 
inertial forces, the boundary layer separates from the body. Depending on 
the length-to-depth ratio the flow field inside the cavity varies (Fig. 2). 
In deeper geometries, the cavity may be bridged by the shear layer which is 
inherently unstable (open cavity.) In shallower geometries, the shear layer 
is deflected inwards with a possible reattachment point at the bottom of the 
cavity (closed cavity). The resonant effects are associated with compressi- 
bility or free surface wave phenomena. The solid boundary may engage in an 
elastic motion coupling with the flow oscillations. 

Longitudinal vortices (Taylor-Gortler like vortices) originate in the 
zone of separation between the primary eddy and the downstream secondary 
eddy, and are convected around the cavity. This separated shear layer is 

continuously pulled into and pushed out of the cavity. When the cavity 
resonates the shear layer is deflected and pumps mass (carrying significant 
momentum) into the cavity. This ingested mass is slowed down by various 
dissipative processes within the cavity. Its presence causes the cavity 
pressure to exceed the freestream pressure. The shear layer is then de- 
flected out of the cavity by the excess pressure and mass is pumped out of 
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the cavity with low momentum. This oscillatory process extracts additional 
freestream momentum during the cycle. The combination of pressure drag (due 
to the cavity walls) and momentum drag (due to the transient eddy cycle) can 
represent a violent oscillating force on the body. Furthermore, this flow- 
field becomes dominated by regions of separated flow within a mixed locally 
subsonic and supersonic flow. 

The shallow cavity flow poses an added complexity to the overall flow 

picture when the shear layer is deflected enough to touch the ceiling of the 

cavity (Fig. 3). The inward bound flow is turned by a series of Mach waves 

which in turn coalesce into a shock. An adverse pressure gradient is gener- 
ated by such compression waves, which causes a secondary separation from the 
ceiling of the cavity. Due to the oscillatory behavior of the flowfield the 
points of separation and reattachment are periodic in time. 

Mathematical Statement 

As it is evident from the results of the experimental studies mentioned 
above, the flow field includes expansions, shocks, separations, recircula- 
tions, and reattachments. The computational simulation of this fluid-dynam- 
ical ly rather complex flow field requires the inclusion of viscous terms, 
modelling of turbulence, regionally clustered gridding, and a shock treat- 
ment mechanism. Fully conservative Navier-Stokes equations are written in 
generalized-curvilinear coordinates to enhance the flow adaptiveness. The 
structure of the turbulence is modelled by the modified k-e model (ref. 3) 
for recirculating flows. Shocks are captured using the divergence of velo- 
city test and treated by the artificial dissipative terms. 

The boundary conditions imposed at the solid walls are the conventional 
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impermeability and no-slip conditions. The wall temperature is assumed to 
be the stagnation temperature (adiabatic wall). The pressure gradient nor- 
mal to the surface is assumed to be zero, however, the gradient parallel to 
the surface is computed from the momentum equation. Downstream and outer 
boundaries are computed from the computation zone values by first-degree 
extrapolation (extrapolation of property gradients). The upstream boundary 
is either assumed to have property profiles computed by a boundary layer 
code, or assumed to be at free stream conditions, thence, a leading edge 
shock. To initialize the computation, the variables all over the computa- 
tion zone are given the values of the upstream boundary (to save computation 
time), with the exception of the points inside the cavity, where the flow is 
assumed stagnant initially. 

General form of the fluid flow equations: 


_^ + V • = 0 (1) 

at 


at 


V 



p 


( 2 ) 


— + V • (E ^ - a V + q) = ? • V (3) 

at 


Perfect Gas Law: 

p = pRT = (y - 1) (E - 1 p f) 

2 


(4) 
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where 


, q=-k (VT), 


(5) 


Newtonian Fluid Stress-Strain Relations: 


a=-pI+T , ?*XV.^I+p[7^+ (V^)^] 


( 6 ) 


The governing fluid motion equations are derived featuring the following: 

- strong conservation form 

- transformed into generalized curvilinear coordinate system (C,n) 

- k-e turbulence modeling or two-layer algebraic turbulence model 

- two-dimensional 

- artificial dissipative terms proportional to the second gradient 

of pressure and temperature, and controlled with blending constants 

- molecular viscosity from Sutherland law 


3Ui 9Fi 3Gi 

— + _ + _ = H 

at ac 3n 


(7) 


where 


Ui = U/J, 

Fi = C(F + Si) 5, + {G + S 2 ) 5y]/J 

(8) 

Si = [F + 

Si) Ox + (G + S 2 ) 'lyl/J 

(9) 
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and 
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( 15 ) 



p = !i [2 (i!; * nf. 1 (in * ii)^] . 1 k (IH . 11) (16) 

P 9x 3y dy 3x 3 3x 3y 3 3x 3y 


u^= p Cy k^/e , = IE _t 


Y-1 Pr4 


(17) 


m 


3/2 


K = II _ 
Y-1 Pr, 


m 


(18) 


P = + P. , K = K + K, 

m t m t 


(19) 


Si • (ax)“ (H) (|u| + a) [2LE. + 5lil] 


3x 




( 20 ) 


Sa-CMV^)(|v|*a) [!l il.Illj 

3y 4P 3^2 4T,^ 


( 21 ) 


The Jacobian an(j the metrics are written as follows and obtained 
nunerical ly from the grid generated for each node: 
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( 22 ) 


j = ^ (^ » ^ ^ 

9(x, y) 


X "y - «y "x ” (’<1 V 


-1 


and 


5x = J(yn) » "y = - 

(23) 

nx = - J (y^) . % “ 

(24) 


Turbulence Models : 

Since the k - e turbulence model is devised for high Reynolds nunber 
number flows, a wall function (Ref. 19) is adopted and used for the very low 
Reynolds numbers regions near the walls (Fig. 4). The procedure of using 
this wall function may be summed up as follows: 


1. velocity profile from previous time step, u=f(y), v=f(y) 


2. k(y) = (. 


8U 


F. N- an 


3- y, • 


(Re )» 




9 “ 20 


4- P Vi / [In E* ReJ , E* = e’^'^'^^v / 


Re. 


,8u 3v^ 

^w+2 ~ ^ (— + — ) 

• 3y 3x w+2 
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5. Assume linear variation of t, i.e. 


^ ^w+2 ■ ^ ^ 

= T„. + - y^) 


‘ w+1 w 


V2 


6. y 


^w+l 


w+1 


7. e 


w+1 w+1 


3u av 

(- +-) 
ay ax w+1 

3/2 
k 


(“m * “t) w+1 


‘Vi - 


8. k . = (p C e 1 

w+1 '■ u w+1-^ 


1/2 


Assume within the viscous sublayer: 


u' = a(t) y 
V = b(t) y 


k = 


2 


near the wall 


The derived equations above are subject to the following initial and 
boundary conditions: 
initial conditions: 

" a given boundary layer profile, ^ = ^(y) 
with ^(y > ff) = T 
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- P = P„. P = Po, 


Compute E and T profiles from V(y) 




1 ' 

(^yR/y-I) 


E = p [-L T 
7-1 




- k, e; 


Dk 



e 

oo 


de 

+ 2 f2 

“ dx 



QO 


where 


f2 = 1 - 0.3 exp[-pk2/ye)2] 


upstream boundary: 

same as initial conditions 
solid wall : 


= 0 


no-slip, impermeable 


(13_) a 0 adiabatic wall, where n is the normal 
9n wall direction to the wall 
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3t 3t 

= + . p 

3X 3X 3y 

Compute 



downstream and outer boundary: 


- first degree extrapolation from the computation zone, i.e. 
gradient of a primitive variable at the boundary is equal to the 
gradient of that variable at the neighboring internal point. 

In view of the upcoming three dimensional simulations, where computer 
storage and CPU time may become limiting constraints, a two-layer algebraic 
turbulence model of Baldwin and Lomax (Ref. 20) is also built in as an 
option. When this option is turned on the k-e equations are 'blanked 
out' . 

Turbulence of the wall shear flow is called the inner turbulence, and 
it is modeled using the Prandtl mixing length theory with the improvement of 
the Van Driest formulation. Using the nomenclature of Baldwin and Lomax, 
inner eddy viscosity coefficient is calculated as 

(»t)lnner-»«' I“1 (25) 

where i s 0.4y [1 - exp (-y /26)] (26) 


12 



and 


(27) 


I 

/ = (P„ y/y^ 

where y is the local normal distance from a solid surface. 

Turbulence of the outer region is based on the Clauser formulation 
using the Klebanoff approximation to the Gaussian error function that gives 
the intermittency factor. The outer eddy viscosity coeficient is calculated 
as 

^^t^outer “ 0*0269 p I^k-i0b (28) 

where the value of Klebanoff intermittency factor, is given by 

fk1eb • * 5.5 (0.3 y/v„,,)‘j-‘ (29) 

and 


F(y) = y |t*)l [1 - exp (- y"^/26)] (30) 

where F(y) is maximum, F , at y = Y . F is the minimum of 

rndx fnsix wdxs 

(^max ^max) O'- (“•25 w where 

= {max [li/p] - min [^/p]} (31) 

at a fixed tangential location. The crossover point, where the switch is 
made from (i^^^outer ^'^t^inner, is the minimuTi value of y where 

("t’outer “ (“t’lmef ^ (“t’outer 

becomes less than (y.).^ 

t 1 nner 
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3. COMPUTATIONAL METHODS 


Explicit- Implicit MacCormack Scheme ; (EMAC): 

Since the mainstream flow is supersonic, the set of partial 
differential equations (7) are predominantly hyperbolic. Therefore, the 
explicit-implicit scheme of MacCormack (Ref. 21), which is developed for 
hyperbolic equations, is one of the methods adopted. The method is second- 
order accurate both in time and space which is achieved through using 
predictor and corrector steps. Each step is consisted of an explicit 
substep and an implicit step. The reasons for this substepping are: (1) 

reduce tridiagonal flux Jacobian matrices to bidiagonal matrices and apply 
direct solution, (2) when time step At meets the stability criterion 
(CFL), the implicit substeps which require most of the computation time are 
automatically bypassed, thence an explicit scheme. 

|Aj and |B| are matrices with positive eigenvalues, related to the 
Jacobi ans 

A = (aF/aO) and B = 36/3U) . 

(a^/a?) and (A^/An) are one-sided forward differences. 

(a /ac) and (a /An) are one-sided backward differences. 

The Jacobi ans A and B are related to the Jacobian A = 3F/3U and 
B = 3G/3U by 


A. 

B = An + Bn 
X y 


(32) 
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1. Explicit Predictor Step: 


A U.. = 
IJ 


- At 


" A-^ F." 


L A5 


a"^ G 


An 


.11 . -1 
1 J 


at every point inside 
the calculation 
domai n 


(33) 


where > (U,/*' - | 


2. Implicit Predictor Step: 


Solve 

'l - At l*i"’ 


I- At a"^|B| 


L J 


An 


= AU." 


ij 


ij 


(34) 


for 6 




n+1 


This is done in two steps. 


Let 


★ 



I - At 




(35) 


This results in two equations, namely, 


a) 


I 


At 


^ n 
A 1A|" 

aF~ 




(36) 
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The above two equations can be solved by splitting |A| to 
$5 Dj Sj-1 as. 


I - at 


+ I 

A |A| 
AS 


★ 


. au,j 


(38) 


Now this equation may be simplified by algebraic manipulation as. 


I - At 




But |A|" . Sj-* 

«0,j 

= AU. . + At A . 1 . 6U- 1 • 

1J ' H+1,0 


IJ 


ii+l,j '“^ij 


^r\ 

(39) 

* 


(40) 


(41) 


n 




Let W= AU,.. . At IAI., 1 ^. 6U^^i^j.then 


6U,.* = (I+AtD^)-l.yl W 


(42) 



This equation can be easily solved since and are known and 

* 

the inversion of the diagonal matrix (l+AtD») is routine. After all 6U.. 

A 1 J 

inside the computation domain are determined the step (b) is carried out in 
the ScBTie manner. Corrector steps are analogous to predictor steps. 

The details of methodology are included in Appendix. 

Line Gauss- Seidel Upwind Relaxation Scheme (IIL6SR) ; 

In the MacCormack scheme explained above, dissipative terms are added 
for nunerical stablity. Often the optimal determination of these terms is 
not straightforward. Furthermore, such terms cause the “smearing" of the 
solution at large-gradient zones and may cause drainage of the conserved 
properties. 

Recently, there has been considerable progress in the development of 
upwind methods, which recognize the hyperbolic nature of the time-dependent 
inviscid equations in the construction of naturally dissipative schemes. 

The improvement in physical treatment comes generally at the expense of 
increased computational work in comparison to central difference approxima- 
tions. However, improved algorithms tailored to the properties of the up- 
wind discretization are being developed. The diagonally dominant properties 
of coefficient matrices arising from such discretizations allow efficient 
relaxation procedures to be developed which can increase the overall 
convergence rate, and thus offset the increased computational work per time 
step (Refs. 22 and 23). 

An efficient relaxation algorithm for the Navier- Stokes equations is 
obtained by using upwind differencing for the convective and pressure terms 
and central differencing for the viscous shear and heat flux terms. The 
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upwind differencing in the present work is implemented using the flux split- 
ting method developed by Van Leer including third-order accurate spacial 
differencing, although the techniques described could be applied to most 
upwind difference methods. 

Rewriting Eq. 7 with slightly different nomenclature is convenient at 
this point. 

+ i£ + iH + ^ + n S)] (43) 

3t 3n 3? 

+ + 5yS)]} 

Q = Q/J 

G = (n^G + (44) 

H = (c G + c H)/J 

^ J 



“p ~ 


~PU ^ 


“pv — 

Q = 

PU 

G = 

PU^+P 

H = 

PUV 


PV 


PUV 


pv^+p 


J- _ 


_u(E+p)_ 


_v(E+pj_ 


where n is the coordinate along the body and ? is the coordinate normal 
to the body. 

The viscous terms on the right-hand side are given by 


~o - 


no - 

T 


T 

XX 


xy 

T„ 

s = 

T 

xy 


yy 

3 _ 


J- _ 


(46) 
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( 47 ) 


’xx ’ ("x \ * 'x^l * * Vc’ 

tyy = (»+2») ("y * ?y \) * »(%“„ * C,U^) (47a) 

'xy “ "‘"yV S”! * * 'x 't> 

“ '^^xx ''’ '^^xy 

+ wPr-i(Y-l)-i (\3j^a2 + C^3^a2) (47c) 


$4 = u"^xy ■'■ ''^yy 

+ uPr'MY-D'Uny 9^ a2 + 9^a2) 


(47d) 


An implicit, upwind, finite-volume scheme described in Refs. 22 and 23 is 
used to solve (1). Application of approximate factorization allows the 
system of equations to be solved in two sweeps; 



+ 3 -A+ + 3 +A- 
n n 

- Re" 3^j" N]aQ* = -RHS 


[— + a; + 3. "^B- 

1 AQ* 

- Re 3 J M]AQ = ^ 

^ jAt 

-+ + - -+ +- 

RHS =3 G +3 G +3 H +3 H 

".1 Jl ^ ' 

-Re {3^[J (n^R + n^S)] 




(48) 


(49) 


(50) 
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The (+) and (-) superscripts indicate flux split quantities, ac- 
cording to the flux vector splittings of Van Leer. A and B are lineariza- 
tions of the G and H fluxes, and M and N arise from the linearizations of 
viscous terms in the c " direction and n - direction, respectively. All 
viscous terms are centrally differenced, and implicit cross-derivative terms 
are neglected in the formulation. 

A comparative study of these two schemes is to be reported in Ref. 24. 

4. RESULTS AND DISCUSSION 

Two Navier- Stokes codes are developed as the tools of analysis. 

Both codes are second order accurate in time and in space. Code A is the 
explicit and time unsplit MacCormack scheme, which proves its use from 
straightforward vectorization point of view (Ref. 21). Code B is an impli- 
cit fully-upwind 1 ine-Gauss-Seidel relaxation scheme that necessitates flux- 
vector splitting (Ref. 22). For a class of unsteady flows, the explicit 
scheme may be more computer-time efficient. However, despite its parallel- 
ization complications due to recursiveness, the implicit scheme is much 
faster converging for steady flows, thence often preferred. 

In an attempt to demonstrate the integrity of these schenes for various 
flow dynamics phenomena and perform comparisons, a variety of high-Reynolds- 
number flow cases are analyzed. These include flow over a flat over a flat 
plate with a leading edge (Fig. 5), flow over a sharp concave corner (Fig. 
6), and flow over a rearward facing step (Fig. 7). Among the phenomena 
observed through simulation are, initiation and growth of a boundary layer, 
shock-boundary layer, shock-boundary layer interaction, expansion around a 
ninety-degree corner and separation, and detached shock. 

The flow simulations are performed using CDC VPS-32 vector processor 
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which has a virtual memory architecture. This manory-to-memory machine has 
short stops which does not use vector registers. For this SIMD (single- 
instruction-multiple-data) processing, the user has to do the memory manage- 
ment and prescribe the required resources, such as, the working set sizes, 
to the operating system. Page memories used by the code are not necessarily 
adjacent, therefore, vector lengths larger than a page (64K-words) are 
avoided. Thirty-two bit (half-precision) arithmetic option and very effi- 
cient scatter-gather operations are available, therefore, they are widely 
used in this code development. The vectors shorter in length than the 
breakeven length are retained in scalar DO loops. The array indeces are 
arranged in decreasing lengths and the nested DO loops containing these 
indeces are ordered likewise, so that, the innermost loop corresponds to the 
longest dimension. Avoidable uses of vector logical unit are eliminated. 

The implicit scheme poses the recurrence problem from parallelization 
standpoint. Most of the operations of the inversion, however, occurs in the 
LU decomposition for four by four block line inversions, and they can be 
computed simultaneously for all the lines before beginning the relaxation 
sweeps. The Gauss-Seidel line inversions constitute less than thirty per- 
cent of total computations as indicated by the histograms (approximately 23% 
in turbulent and 26% in laminar calculations). Therefore, a further im- 
provement in computation times would have been achieved by employing a par- 
allel algorithm. The scalar computational time on VPS-32 for the above 
algorithm is approximately 800 microseconds per gridpoint per time step. 

This vectorized upwind scheme achieves a speed up to approximately 17 micro- 
seconds per gridpoint per time step for laninar Navier Stokes calculations 
and 21 microseconds for turbulent Navier Stokes calculations. Another 21% 
time reduction can be obtained if thin-layer Navier Stokes equations are 
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used. The scaler computation time for the MacCormack code, 240 microseconds 
per grid point per time step, is reduced to 40 microseconds through vectori- 
zation. Better speed-ups are possible for larger mesh sizes. 

Further discussion on computational methodology comparisons are dis- 
cussed in Ref. 24. This reference also includes analyses and simulations 
for flow over an airfoil and flow over an ogive-nosed cylinder at angles of 
attack. Such cases are of interest to the second part of this investiga- 
tion, therefore, not included in this report. 

Two-dimensional simulations of cavity flowfields are demonstrated in 
Figures 8 through 25. These results certainly do not reflect the third- 
dimension effects of the side walls. The box-cavity flowfield is synmetri- 
cal with respect to its centerplane which is at z = width/z (Fig. 1). Two- 
dimensional simulations reflect the flowfield at this plane of synmetry. 
Example cases included are two open cavities (length-to-depth ratios are 3 
and 6), a transitional cavity (length-to-depth ratio is 12), and two closed 
cavities (length-to-depth ratios are 16 and 24). Figures corresponding to 
these cases are 11 through 13, 14 through 16, 17 through 19, 20 through 22, 
and 23 through 25, respectively. All of these figures display the final 
property values at 0.1 seconds after the computations are intiated. These 
results are also reported in References 16 through 18, which are authored or 
co-authored by the principal investigator of this research. 

5. CONCLUDING REMARKS 

An explicit-implicit and an implicit two-dimensional Navier-Stokes 
code along with various grid generation capabilities are developed. A 
series of classical benchmark cases are simulated using these codes. 
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References 16, 17, 18, and 24. Current follow-up efforts may be sunmed up 
as; developing explicit as well as implicit three dimensional codes, three- 
dimensional cavity flowfield studies and simulations, two-dimensional and 
three-dimensional store flowfield studies and simulations. 


23 


REFERENCES 


1. Stallings, R. L., Jr., "Store Separation from Cavities at Supersonic 

Flight Speeds," Journal of Spacecraft and Rockets, Vol. 20, No. 2, 
pp. 129-132, 198T ^ 

2. Blair, A. B. and Stallings, R. L., "Supersonic Axial-Force Characteris- 

tics of a Rectangular-Box Cavity with Various Length-to-Depth 
Ratios in a Flat Plate," NASA-TM-87659, April 1986. 

3. Rossiter, J. E., "Wind Tunnel Experiment on the Flow Over Rectangular 

Cavities at Subsonic and Transonic Speeds," R&M No. 3438, British 
A.R.C., October 1964. 

4. Charwat, A. F., Rous, J. N., Dewey, F. C., and Hitz, J. A., "An 

Investigation of Separated Flows," J. of Aerospace Sciences, vol. 

28 Part I, June 1961, Part II, July 1961. 

5. McDearmon, R. W., "Investigation of the Flow in a Rectangular Cavity in 

a Flat Plate at a Mach Number of 3.55," NASA TN 0-523, September 
1960. 

6. McGregor, 0. W., and White, R. A., "Drag of Rectangular Cavities in 

Supersonic and Transonic Flow Including the Effects of Cavity 
Resonance, "AIAA Journal, vol. 8, November 1970. 

7. Sinha, S. N., Gupta, A. K., Oberai, M. M., "Laminar Separating Flow 

Over Backsteps and Cavities," AIAA Journal, Vol. 20, No. 3, 1982. 

8. Antonov, A. N. and Shalaev, S. P., "Experimental Investigation of 

Nonstationary Flow in Cavities in a Supersonic Flow," Fluid 
Dynamics-Soviet Research, vol. 14, no. 5, 1980. 

9. Varfolomeyev, I. M., Gortysher, Y. F. and Shekukin, V. K., "Heat 

Transfer and Dynamics of Supersonic Airflow Past Cavities," Heat 
Transfer-Soviet Research, vol. 13, 1981. 

10. Clark, R. L., Kaufman, L. G., Maciulaitis, A.>, "Aeroacoustic 

Measurements for Mach 0.6 to 3.0 Flows Past Rectangular Cavities," 
AIAA Paper No. 80-0036, 1980. 

11. Borland, C. J., "Numerical Prediction of the Unsteady Flowfield in an 

Open Cavity," AIAA Paper No. 79-0136, 1979. 

12. Hankey, W. L.; and Shang, J. S., "Numerical Solution to Pressure 

Oscillations in an Open Cavity," AIAA Paper No. 79-0136, 1979. 

13. Gatski, T. B. and Grosch, C. E., "Embedded Cavity Drag in Steady 

Laminar Flow," AIAA Journal, Vol. 23, No. 7, 1985. 

14. Gorski, J. J. and Chakravarthy, S. R., "Calculation of Three 

Dimensional Cavity Flowfields," AIAA Paper No. 87-0117, 1987 (to be 
presented) . 


24 



15. dm, D., "Navi er- Stokes Simulation for Flow Past an Open Cavity," AIAA 

Paper No. 86-2628, 1986. 

16. Baysal, 0., Stallings, R. L., Srinivasan, S., "Unsteady Flow in a 

Cavity at Supersonic Speeds and its Navier- Stokes Prediction with a 
Vectorized Algorithm," AIAA Paper No. 86-0594, 1986. 

17. Baysal, 0., Stallings, R. L., and Srinivasan, S., "Cavity Flowfield 

Calculations," Virginia Journal of Science, vol. 37, no. 2, p. 39, 
Summer 1986. 

18. Baysal, 0., Stallings, R. L., "Computational and Experimental 

Investigation of Cavity Flowfields," AIAA Paper No. 87-0114, 

January 1987 (to be presented) . 

19. Sindir, M. M. S., "A Numerical Study of Turbulent Flows in Backward 

Facing Step Geometries," Ph.D. Thesis, University of California, 
Davis, 1982. 

20. Baldwin, B. S. and Lomax, H., "Thin Layer Approximations and Algebraic 

Model for Separated Turbulent Flows," AIAA Paper No. 78-257, 1978. 

21. MacCormack, R. W., "Current Status of Numerical Solutions of the 

Navier- Stokes Equations," AIAA Paper 85-0032, 1985. 

22. Thomas, J. L. and Walters, R. W., Upwind Relaxation Algorithms for the 

Navier- Stokes Equations. AIAA Paper No. 85-1501 CP, July 1985. 

23. Rumsey, C. L., Thomas, J. L., Warren, 6. P., and Liu, G. C.,, "Upwind 

Navier-Stokes Solutions for Separated Periodic Flows," AIAA Paper 
No. 86-0247, 1986. 

24. Baysal, 0., "Flow Simulations Using Upwind Line-Gauss-Seidel and 

MacCormack Schemes On a Vector Processor," 1987 ASME Fluids 
Engineering Conference, Cincinnati, Ohio, June 1987 (to be 
presented) . 

25. Carter, J. E., "Numerical Solutions of the Navier-Stokes Equations for 

the Supersonic Laminar Flow Over a Two-Dimensional Compression 
Corner," NASA-TR-R385, July 1972. 


25 





L 


C^) 


Figure 2. Hypothetical structure of flowfield variation due to depth-to 
length ratio: (a) Closed or shallow cavity, (b) open or deep 
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Figure 5. continued 
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Figure 8. Sample grid meshes for cavity flowfield calculations: (a-b) body- 
fitted meshes, (c-d) patched rectangular meshes. 
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Figure 9. Continued 
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Figure 10. Experimentally determined pressure distributions of cavity 
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Figure 11. Supersonic flow over L/D = 3.0 cavity: Pressure coefficient 

distribution at, (a) front flat plate, Spp]^, (b) forward face, 
Sep, (c) ceiling, Sq, (d) rearward face, Snp, (e) rear flat 
plate, Spp 2 » (f) skin friction coefficient, (g) displacement 
thickness, (h) momemtum thickness distributions. 
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Figure 13. Supersonic flow over L/D = 3.0 cavity: (a) Pressure contours, 
(b) density contours, (c) Mach hunber contours. 
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Figure 15. Supersonic flow over L/D = 6.0 cavity; (a) streamline contours, 
(b) velocity vectors, (c-e) blow ups, of corner flows. 
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Figure 16. continued 
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Figure 16. Supersonic flow over L/D = 6.0 cavity: (a) Mach number contours, 
(b) Pressure contours, (c) Density contours, and (d) total 
pressure loss contours. 
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Figure 17. continued 
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Supersonic flow over L/0 = 12.0 cavity: (a) pressure coefficient 
distribution over entire wall, (b) skin friction coefficient 
distribution over entire wall, (c-g) pressure coefficient 
distribution on front flat plate (Sppj^), forward face (Sep), 
ceiling (S_), rearward face (S|^p), rear flat plate (Spp 2 ;, 
respectively (see Fig. 10 for comparison). 
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Figure 18 . continued 
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Figure 18. Supersonic flow over L/D = 12.0 cavity; (a) velocity vectors at 
an earlier stage of flow development, (b-c) streamlines and 

velocity vectors of developed flowfield, (d-e) Blow ups of 
corner flows. 
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Figure 19. Supersonic flow over L/0 = 12.0 cavity: (a) pressure contours, 
(D) density contours, (c) Mach number contours. 
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Figure 21. Supersonic flow over L/D = 16.0 cavity: (a) streamline contours 
(b) velocity vectors of complete flowfield, (c-d) Blow ups of 
corner flows. 
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Figure 22. Supersonic flow over L/D = 16.0 cavity: (a) pressure contours 
(b) density contours, (c) Mach nunber contours. 
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Figure 23. continued 
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Figure 23. Supersonic flow over L/D = 24.0 cavity; (a) pressure coefficient 
distribution over entire wall, (b) skin friction coefficient 
distribution over entire wall, (c-g) pressure coefficient 
distribution on front flat plate (Sppj^), forward face (Sep), 
ceiling (S..), rearward face (Snp), rear flat plate (Spp 2 )> 
respectively. (See Fig. 10 for comparison). 
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Figure 25. Supersonic flow over L/D =* 24.0 cavity: pressure contours, (b) 
density contours, (c) Mach nunber contours. 
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APPENDIX 

DETAILS OF MACCORMACK SCHEME ADAPTATION 
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The integration scheme can be simplified by di agonal ising A and B 



Once the eigenvalues of A and 8 are known, it is possible to express 
and B in the form 



where and Ag are diagonal matrices consisting of the eigenvalues of 

A, A^^j^ •••^A,4 ^8,1’ •••^8 4’ '^®spectively. The vectors 

and are constructed using the eigenvalues of A and 8, respectively, 
as columns. 

where 
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Dg = |Ag| + Ag I. 


where 
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B 



Ati 

At 


V = max [y, (a+2u) , k] 


The expressions for 0^ and Dg assume that viscous effects are modeled in 
the implicit part of the scheme by the addition of the terms A^, Ag which 
include viscosity through the coefficient, v. 
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